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Ionization of highly charged relativistic ions by neutral atoms and ions is considered. Numerical 
results of recently developed computer codes based on the relativistic Born and the equivalent-photon 
approximations are presented. The ionization of the outer shells dominate. For the outer projectile 
electron shells, which give the main contribution to the process, the non-relativistic Schrodinger wave 
functions can be used. The formulae for the non-relativistic reduction of the Dirac matrix-elements 
are obtained for ionization of electrons with arbitrary quantum numbers n and I. 



PACS numbers: 34., 34.10.+X, 34.50.-s 

I. INTRODUCTION 

Projectile ionization, also referred to as electron loss or 
stripping, is an important charge-changing process play- 
ing a critical role in many applications such as heavy-ion 
driven inertial fusion, ion-beam lifetimes in accelerators, 
medical science, material technology and others (see. e.g., 

H-B). 

In order to calculate the non-relativistic ionization 
cross sections for ion-atom collisions a computer code 
called LOSS was developed using the Born approxima- 
tion [3] with account for the screening as well as antis- 
creening effects. Non-relativistic Schrodinger wave func- 
tions are used to describe the electronic structure. This is 
a good approximation for the outer shell electrons, which 
play the dominant role. 

Relativistic ionization has been pioneered at Stanford 
by Anholt and collaborators [5[ and it is reviewed in sev- 
eral books and review articles (see e.g., @- @). Numer- 
ical calculations of relativistic ionization cross sections 
are presented in the literature mainly for ionization of 
H- and He-like ions from the ground Is state (see, e.g., 

0-0)- 

Recently, two new computer codes LOSS-R (Relativis- 
tic Loss) and HERION (High Energy Relativistic ION- 
ization) have been developed for calculation of the rela- 
tivistic ionization cross sections. The LOSS-R code (l6j 
was created on the basis of the non-relativistic LOSS 
code [|| using the relativistic Born approximation in the 
momentum-transfer representation without magnetic in- 



teractions. The HERION code [ijj uses the dipole and 
impulse approximations with relativistic wave functions 
but it also neglects the magnetic interactions. 

It is the aim of the present work to provide a theo- 
retical framework for ionization in relativistic ion-atom 
collisions with possible account for the magnetic interac- 
tions for an arbitrary many-electron ions. The formulae 
obtained can be used for projectile ions such as U 28+ col- 
liding with the rest-gas atoms and molecules at energies 
up to a few GeV/u. Such heavy many-electron ions are of 
practical implications for the international FAIR project 
at GSI Darmstadt [H[. 

In section 2 the main formulae for relativistic ion- 
ization are recapitulated, also in order to establish the 
notation. Numerical results due to the LOSS-R and 
HERION codes are given in section 3. In section 4 the 
non-relativistic reduction of the Dirac matrix-elements is 
done. 



II. FIRST ORDER PERTURBATION THEORY 
OF RELATIVISTIC IONIZATION 

An accurate procedure to include the relativistic ef- 
fects in ionization cross sections using the plane wave 
Born approximation is given in [||. There, the cross sec- 
tion summed over the magnetic quantum numbers and 
presented in the momentum-transfer ^-representation is 
given in the form: 
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where a = e 2 /he denotes the fine-structure constant, Zt 
the effective charge of the target atom, q — lu/v, and uj 
is the ionization energy. The ion velocity is given by the 
relativistic factor (3 — v/c. The matrix-elements F and G 
are evaluated in a x, y, z coordinate system in which the 
z-axis lies along the momentum transfer vector q. The 
x-axis lies in the plane formed by the vectors q and the 
projectile velocity v. They contribute incoherently to the 
cross section due to different selection rules for the final 
magnetic substates. The matrix elements F and G are 
given by Q 



F(q)=<f\e i " z \i>, 



G(q) ^< f\c 
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where \i > and \f > denote the wave functions in the 
initial and final states and a x the x-component of the 
Dirac matrix vector a. This is a useful splitting: the F- 
term tends to a constant for 7 — > 00 whereas the G-term 
increases as Iwy due to a contribution of the 'equivalent' 
photons. The matrix element F is the main term in the 
non-relativistic ionization, e.g. in the formulae of the 
LOSS code Q whereas the matrix element G appears 
for relativistic ionization. The G-term can be calculated 
numerically if the Dirac wave functions with the large 
and small components are known. 

The equivalent photon approximation can be obtained 
from eq. (1): Due to the photon pole, the second term 
will dominate at high energies. For small g-values the 
quantity G(q) is directly related to the dipole matrix el- 
ement: 



G (°) = -jjT < Mi >s -jp D f* 



(4) 



The integral over q in eq. (1) has a dominant contri- 
bution at q ~ qo. In the integration over q, we may 
assume that G(q) = G(qo). Then the integral diverges 
logarithmically at large q. One can introduce a suitable 
cutoff q max which corresponds to a cutoff impact param- 
eter b m in in the impact-parameter space. In turn, the 
dipole moment determines the photoionization cross sec- 
tion (see e.g., eq. (69.2)): 
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So the ion-atom ionization cross section can be expressed 
via a photoionization cross section and the equivalent- 
photon number. 

The PWBA (Plane- Wave Born Approximation) for- 
malism is also used in the calculation of bound-free pair 
production in antiproton-nucleus collisions (antihydrogen 
production) [2(| and in heavy ion collisions [2l]]. Pair 
production may be viewed as ionization of the negative 
energy Dirac sea. In the pair production process one has 
to use Dirac wave functions, rather than the Schrodinger 
wave functions. 

In (l7j the dipole approximation was realized as a part 
of the HERION code where the dipole part of relativistic 
ionization cross section of the projectile ion in collision 
with a neutral atom is expressed via photoionization cross 
section a (to) (see pj): 
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where io m i n denotes the ionization threshold and n{u) 
the equivalent-photon number 
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where K n (x) is the McDonald function. The parameter 
bmin is defined by the size of the projectile-ion shell nl 
from which ionization occurs and is given with a good 
accuracy by 



where ao = 0.53 x 10~ 8 cm is the Bohr radius, I n i denotes 
the binding energy of the nl shell and Ry — 13.606 eV is 
the Rydberg energy unit. 



(8) 
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FIG. 1: Calculated ionization cross sections of 
U 28+ (ls 2 ...4/ 14 5s 2 5p 6 ) ions by protons: the LOSS-R - 
code [la ]. Contribution of ionization from different subshells 
of U 2 are shown together with the total cross section. 



FIG. 2: Relativistic dipole, non-dipole and total ionization 
cross sections, i.e. summed over all nl shells, of U 28+ by pro- 
ton impact calculated by the HERION and LOSS-R codes (in- 
dicated): dashed curves - the non-dipole parts, dotted curves 
- the dipole parts, and solid curves - the total cross sections. 



III. NUMERICAL CALCULATIONS 



Calculated partial (on projectile electron subshells nl) 
and total ionization cross sections of U 28+ ions by pro- 
tons at energies E = 1 - 100 GeV/u are shown in Fig. Q] 
Calculations were performed by the LOSS-R code includ- 
ing dipole and non-dipole parts of the ionization cross 
sections and using the non-relativistic wave functions. 
As seen from the figure, the main contribution to the 
total cross section is given by electron ionization from 
3c? 10 , 4p 6 , 4d 10 ,..., shells which can be described by the 
non-relativistic wave functions 

A comparison of relativistic ionization cross sections 
of U 28+ ions by proton impact calculated by two dif- 
ferent codes, LOSS-R and HERION, is shown in Fig. EU 
where the contributions from dipole-, non-dipole parts 
are shown together with the total (dipole + non-dipole) 
cross sections. In the LOSS-R code, the dipole part of the 
cross sections corresponds to electron transitions into the 
continuum with the orbital angular momentum quantum 
numbers A = I ± 1 where I is the angular momentum of 
the projectile-electron shell. As seen from the figure, the 
non-dipole parts calculated by the codes agree within 30 
% meanwhile the dipole parts and the total cross sec- 
tions agree only within a factor of 2. Most probably, this 
discrepancy is related to two reasons: first, the use of dif- 
ferent wave functions in the codes, i.e., relativistic wave 
functions in HERION code and and non-relativistic ones 
in the LOSS-R code. And second, in the HERION code 
dipole ionization cross sections are calculated on the ba- 
sis of photoionization cross sections; the latter are much 
higher than those calculated in the usual non-relativistic 



approximation because of the presence of the so-called gi- 
ant resonances occurring due to configuration interaction 
of the final states (see [22| for details). 

We note that the dipole part of ionization cross section 
involves about 60-70 % of the total cross section while 
the relativistic non-dipole part has a weak dependence 
on energy and, therefore, the dipole part has practically 
the same shape as the total cross section. 

Relativistic ionization cross sections of Au 78+ (ls) ions 
in collisions with carbon are displayed in Fig. [3] where ex- 
perimental data are given by symbols and theoretical cal- 
culations by curves. The non-relativistic result (LOSS) 
shows a Born maximum followed by a decrease of the 
cross section. The relativistic result (LOSS-R) has a lo- 
cal minimum around 1 GeV/u energy and increases log- 
arithmically with energy. The results obtained with the 
HERION code have a little better agreement with ex- 
periment at E > 200 MeV/u than those by the LOSS-R 
code. In the calculations, the effective charge of the tar- 
get Z T = Z 2 + N = 42 was used while in the LOSS-R 
code the Zt value is found in the g-represcntation with 
account of all subshells of the target (see eq. 27 and Rcf. 
[3] for details). The lower curve is the recent result [l5j] 
with relativistic electron description. 
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FIG. 3: Ionization cross sections of H-like Au (Is) ions by 
carbon atoms. Experiment: open circles - 12], 23] solid cir- 
cle - HI]. Theory: LOSS - non-relativistic LOSS code Q, 
LOSS-R - relativistic LOSS-R code QJ], HERION - HERION 
code [13] with effective targ et charge Z%- = Z 2 + N = 42, V - 
calculations by Voitkiv [15( with relativistic electron descrip- 
tion. 



IV. NON-RELATIVISTIC LIMIT OF a-MATRIX 
ELEMENT 

Now we consider the non-relativistic limit of the 
matrix-element G given by eq. (|3"|): 

G(q) =< f\a x exp{iqz)\i > . 

The operator J = a x exp (iqz) is an odd one and con- 
nects the large and small components of the Dirac wave 
function. It gives zero if one only uses the large compo- 
nents. The non-relativistic limit of G(q) is given in [251 ] 
in the form: 

J = ■ 7 ^—(a x exp(iqz)p 1 (a-p) + (a-p)p 1 a x exp(iqz)), (9) 
2m 

where p = e z q. This operator acts only on the large com- 
ponent of the Dirac wave function. We use the relations 
ct x pi = <J X and piOi x — & x , (see [25| , p. 891) to obtain 
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(a x exp (iqz)(d ■ p) + (a ■ p)a x exp (iqz)). (10) 



Using the identity (<? • A) (a ■ B) = (A ■ B) + id ■ (A x B) 
(see e.g. eq. XIII. 83 of (2f|) with A = e x ;B = p and 
with B = e x ; A = p we find 

J = - — (exp (iqz)p x + p x exp (iqz) + spin — flip terms). 
2m 

(11) 

The spin-flip terms are given by 

exp (iqz)(p x o% — (fx o% exp (iqz) (12) 

In the following, we neglect the spin-flip terms because 
they are small in our case, see also J5|). So the non- 
relativistic limit of the odd operator J is 



J = —(exp (iqz)p x + p x exp (iqz). 
2m 



(13) 



This result can be also obtained directly using the 
electromagnetic-current approach in the Schrodinger the- 
ory (see, e.g., (2(|, Chapter 82a). 

The spin-flip terms appear in the Pauli approximation, 
which can be obtained from the Dirac equation. The 
matrix element G(q) becomes 



G(q) = ^<f\eM*qz)^> 



(14) 



Here p x = —i^4z and 



^iqz 



0. 



The derivative 



is easier to calculate in spherical 
d_ 



dz 

coordinates rather than so we rotate the coordinate 



system by 7r/2 around the y-axis. This changes x — > z 
and z — > —x. So we have 

il% d 
G(q) = — < f\exp(-iqx) — \i> . (15) 
mc oz 

Actually we only need the absolute square of the func- 
tion G, summed over spins and integrated over the angle 
of the outgoing electron. This matrix element is calcu- 
lated according to the rules of angular momentum alge- 
bra (see e.g., [23]). The bound state wave function \i > 
is presented in the form \i >— j\ i (r)Yi iTni (f), i.e. with 
separated angular and radial parts. 

The exponential factor is expanded into partial waves 
in the usual way: 



exp (-iqx) = A-K^i +l ji(qr)Y lm (r)Y? m (i 



(16) 



The radial wave function of the bound and continuum 
states satisfy the normalization conditions: 



9u ( r ) dr = 1' 91 f M ~ sin 
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-ln(2fcr) 
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k 2 /2 = e, 



(17) 



where e is the energy of ejected electron and r/ is the scattering phase shift. 
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According to 27] , the action of d/dz on this wave 



Thus, the matrix element is presented as a sum of angular 
and radial parts. We note that the present procedure is 
quite similar to that used for calculating the magnetic 
terms of the pionium breakup [28| . 

We define two types of radial matrix elements, a 'usual' 
one 

/>oo 

RffXlM = / drg lfJX (qr)fu, (19) 
Jo 

and a 'new' one involving the derivative of the initial 



Collecting all the 



G 



x 



As was mentioned above, we need the function |G(<7) 2 | 
integrated over fi^ and summed over spins. The integra- 
tion over makes the sum over If and to/ incoherent, 
due to the orthogonality of the spherical harmonics. We 

I 



function yields: 



(18) 



radial wave function: 

Rt fXH (q(q) = jH drrg lfjx (qr)^. (20) 

The three-dimensional integration is, as usual, decom- 
posed into a radial and angular parts. The angular 
integration can be done using eq. (4.6.3) in [13] and 

In the integration over the three spherical harmonics 
one obtains two types of terms: f2+ and f2„ with 



(21) 



(22) 



assume that the states j = Z± 1/2 are degenerate. Using 
the sum rules, the summation over to^ and mj can be 
easily done: 
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The summation over can be also done using the com- pleteness relation 
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and one obtains the function |G(q))| 2 for given projectile quantum numbers rij, U 
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The radial integrals R B and R d depend on q, l/,X, k 
according to eqs. (19) and (20), and u± is given by 



We note that A = corresponds to the dipole exci- 
tation, and a monopole contribution does not exist for 



J 



the matrix element G but corresponds to the 'equivalent 
photon' contribution with the photon having spin 1. The 
monopole part is included to the the matrix-element F 
(see eqs. (1) and (2)). 

Finally, for relativistic ionization cross section from the 
projectile electronic shell by a heavy target particle one 
can use eq. (1) but replacing the target effective charge 
Zt by the g-dependent charge Zx{q) in the form: 



N T 



Z T (q) = [Z-^F^q? 



N t -Y. F iM 2 ] - Fjj(q) =< j|exp(iqr)|j >, 



(27) 



r 



where Z and Nt denote the nuclear charge and the num- 
ber of electrons of the target [Z = Nt for neutral atoms 
and Z = 1, N T = for protons). The form of eq. [23 
makes it possible to take into account the screening and 
anti-screening effects of all target electrons. 

The function F(q) in eq. (1) is the 'usual' Born matrix 
element which after separating angular and radial parts 
is written in the form: 

™ 2 | = (2A + 1)(2Z / + 1)(^ */ J) 2 |i? s (g)| 2 , (28) 



where the Born integral R (q) is given in eq. JT9 



The function G(q) describes the magnetic interactions 
between projectile ion and target atom and has the form: 



|G(g)| 



2 2 (2A + 1)(2Z / 
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If A k-1 




[kR d (q)+h(h-l)R B (q)} 



(39) 



where the radial integral R d {q) is given by eq. (|20|) . 



V. CONCLUSION 

Relativistic ionization of heavy ions colliding with 
atoms or ions is considered in the plane-wave Born ap- 



proximation (PWBA). The ionization cross sections are 
presented in the momentum-transfer ^-representation as 
a sum of two terms: the 'usual' Born approximation (ma- 
trix element of iqr) and the relativistic term responsible 
for the magnetic interactions between colliding particles 
and expressed through the x-component a x of the Dirac 
matrix vector a. 
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A simple limit of PWBA is the equivalent photon 
approximation (dipole approximation). Numerical cal- 
culations for the dipole and non-dipole parts of rela- 
tivistic cross sections are presented for the p + U 28+ 
and C+Au 78+ collisions using the LOSS-R and HERION 
computer codes. 

A new formula is obtained for the relativistic ionization 
cross section for an arbitrary nl shell of the projectile 
ion in the form of separated angular and radial parts. 
The radial part is expressed via an integral of the initial 
bound state wave function and the derivative of the final 
continuum wave function. In this way, new possibilies to 



investigate ionization processes including many electron 
atoms are opened. 
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